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We discuss a recently developed formalism which describes 
the quantum evolution of a solid-state qubit due to its contin- 
uous measurement. In contrast to the conventional ensemble- 
averaged formalism, it takes into account the measurement 
record and therefore is able to consider individual realizations 
of the measurement process. The formalism provides testable 
experimental predictions and can be used for the analysis of 
a quantum feedback control of solid-state qubits. We also 
discuss generalization of the Bayesian formalism to the con- 
tinuous measurement of entangled qubits. 



I. INTRODUCTION 

Bayesian approach to the problem of continuous quan- 
tum measurement is a relatively new subject in solid- 
state mesoscopics, even though this approach has a long 
history [1,2] as a general quantum framework and is 
rather well developed, for example, in quantum optics 
[3] (for more references, see Ref. [4]) The main problem 
considered in this paper is a very simple question: what 
is the evolution of a quantum two-level system (qubit) 
during the process of its measurement by a solid-state 
detector (Fig. 1)? In spite of the question simplicity, the 
answer is not that trivial. 

The textbook "orthodox" quantum mechanics [5] says 
that the measurement should instantly collapse the qubit 
state, so that after the measurement the qubit state is ei- 
ther |1) or 1 2), depending on the measurement outcome. 
[The measurement basis is obviously defined by the de- 
tector; in particular, it is a charge basis for the exam- 
ples of Figs. 1(a) and 1(b).] Such answer is sufficient 
for typical optical experiments when the measurement 
is instantaneous (a scintillator flash or a photocounter 
click). However, for typical solid-state setups (as well as 
for some more advanced setups in quantum optics [3]) the 
instantaneous collapse is not a sufficient answer. In par- 
ticular, in the examples of Fig. 1 typically the detector 
is weakly coupled to the qubit, so the measurement pro- 
cess can take a significant time and therefore the collapse 
should be considered as a continuous process. The no- 
tion of a continuous evolution due to measurement is well 
accepted in the solid-state community and is usually con- 
sidered within the framework of the Leggett's formalism 
[6,7]. This formalism gives the decoherence-based answer 
to the question posed above. It says that the nondiagonal 
matrix elements of the qubit density matrix (obtained by 
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FIG. 1. (a) General schematic of a continuously measured 
solid-state qubit and two particular realizations of the setup: 
(b) a qubit made of double quantum dot (DQD) measured 
by a quantum point contact (QPC) and (c) a qubit based on 
single-Cooper-pair box (SCPB) measured by a single-electron 
transistor (SET). 



tracing over the detector degrees of freedom) gradually 
decay to zero, while the diagonal matrix elements do not 
evolve (assuming that the qubit docs not oscillate by it- 
self, iJ = 0, where H describes the tunneling between |1) 
and 1 2)). So, after the completed measurement we have 
an incoherent mixture of the states |1) and |2). 

Let us notice that these two answers to our question 
obviously contradict each other and the "orthodox" an- 
swer cannot be obtained as some limiting case of the 
decoherence answer (since decoherence does not lead to 
localization into one definite state). The resolution of the 
apparent contradiction is simple: two approaches con- 
sider different objects. The decoherence approach de- 
scribes the average evolution of the ensemble of qubits, 
while the "orthodox" quantum mechanics is designed to 
treat a single quantum system. This difference also ex- 
plains the inability of the decoherence formalism to take 
the measurement outcome into account. 

Obviously, it is desirable to have a formalism which 
would combine advantages of the two approaches and 
describe the continuous measurement of a single qubit. 
Then the "orthodox" result would be a limiting case for 
very fast (and "strong") measurement, while the deco- 
herence result could be obtained by an ensemble aver- 
aging. The Bayesian formalism [4,8] which is the sub- 
ject of this paper has been developed exactly for that 
purpose (some extensions of the Bayesian formalism will 
be discussed later). Notice that the formalism has been 
also reproduced in a somewhat different language (us- 
ing the terminology of quantum trajectories, quantum 
jumps, and quantum state diffusion) by another group 
[9]. It is important to stress that the Bayesian approach 
is not a phcnomcnological formalism which just correctly 
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describes two previously known cases. It claims the de- 
scription of a real and experimentally verifiable evolution 
of a single qubit in a process of measurement. 

Simply speaking, the Bayesian formalism gives the fol- 
lowing answer to the question posed above (for H = 0). 
During the measurement process the diagonal matrix ele- 
ments of the qubit density matrix evolve according to the 
classical Bayes formula [10,11] which takes into account 
the noisy detector output [I{t) in Fig. 1] and describes a 
gradual qubit localization into one of the states |1) or |2), 
depending on I{t). The evolution of nondiagonal matrix 
elements can be easily calculated using somewhat sur- 
prising result that a good (ideal) detector preserves the 
purity of the qubit state, so that the decoherence is actu- 
ally just a consequence of averaging over different detec- 
tor outcomes I(t) for different members of the ensemble. 
(Nonideal detectors also produce some amount of qubit 
decoherence, which is calculated within the formalism.) 

Notice that in the case of an ideal detector, our re- 
sult can actually be considered as a simple consequence 
of the so-called Quantum Bayes Theorem (we borrow 
this name from the book on quantum noise by Gardiner 
[12], even though it is not a theorem in a mathemati- 
cal sense). However, the application of this "theorem" 
is not always straightforward, so instead of applying it 
as an ansatz, we derive the Bayesian formalism for par- 
ticular measurement setups, starting from the textbook 
quantum mechanics. 

It is difficult to avoid philosophical questions discussing 
a problem related to quantum measurements. In brief, 
philosophy of the Bayesian approach is exactly the phi- 
losophy of the "orthodox" quantum mechanics. A minor 
technical difference is that instead of assuming instanta- 
neous information on measurement result corresponding 
to instantaneous "orthodox" collapse, we consider a more 
realistic case of continuous information flow. 

Finally, let us mention that the problem of solid-state 
qubit evolution due to continuous measurement was re- 
cently a subject of theoretical study by many groups (see, 
e.g. [13 18]). However, most of them assumed ensem- 
ble averaging and so obtained results different from the 
Bayesian results (except the Australian group [9,19,20] 
which also studies single realizations of the measurement 
process). 



II. SIMPLE MODEL 

Even though Bayesian approach is applicable to a 
broad range of measurement setups, let us start with a 
particularly simple setup [Fig. 1(b)] consisting of a double 
quantum dot occupied by a single electron, the position 
of which is measured by a low-transparency Quantum 
Point Contact (QPC) or (which is almost the same) by 
just a tunnel junction [Fig. 2(a)]. Basically following the 
model of Ref. [13] we assume that the detector barrier 
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FIG. 2. (a) Tunnel junction as a detector of the dou- 
ble-quantum-dot qubit. The electron location in the DQD 
affects the detector barrier height. The noisy current I{t) (de- 
tector output) reflects the evolution of the qubit density ma- 
trix pij(t). (b) Idea of the Bayesian formalism derivation via 
Bloch equations. The number n of electrons passed through 
the detector is periodocally collapsed (forced to choose a def- 
inite value) at moments tk- 



height depends on the location of the electron in either 
dot 1 or 2; then the current through the tunnel junction 
(which is the detector output) is sensitive to the electron 
location. 

The Hamiltonian of the system. 



W = Hqb + Ti-DET + Ti-INT, 



(1) 



consists of terms describing the double-dot qubit, the de- 
tector, and their interaction. The qubit Hamiltonian, 

'Hqb = I (4^2 - 4ci) + H (4c2 + 4ci), (2) 

is characterized by the energy asymmetry e between two 
dots and the tunneling strength H (we assume real H 
without loss of generality). The detector and interaction 
Hamiltonians can be written as 

Ti-DET = ^ Eiajai + ^ Erulur + ^ T{alai + ajur), 



Hint = ^ -y {c\ci - 4c2)(4«i + «;«r), 



(3) 



where both T and AT are assumed real and their de- 
pendence on the states in electrodes (/,r) is neglected. 
For simplicity we assume zero temperature (Bayesian 
formalism at finite temperatures has been considered 
in Refs. [4,9,21,22]). If the electron occupies dot 1, 
then the average current through the detector is Ii = 



2 



27r(r + AT/2fpipre'^V/n {V is the voltage across the 
tunnel junction and pi^r are the densities of states in the 
electrodes), while if the measured electron is in the dot 
2, the average current is h = 2tt{T - bJ: j^Y pip^-e^V jU. 
The difference between the currents, 



AI = h- h 



(4) 



determines the detector response to the electron position. 
Because of the detector shot noise, the two states cannot 
be distinguished instantaneously and the signal-to-noise 
ratio (S/N) gradually improves with the increase of the 
measurement duration. The S/N becomes close to unity 
after the "measurement" time 



2(A/)2 



(5) 



where the spectral densities 5i and S2 of the detector 
shot noise for states |1) and |2) are given by the Schottky 
formula, 



5'i,2 = 2e7i,2. 



(6) 



[Actually, Eq. (5) exactly corresponds to S/N=l for 
= ^2, while for S/N^l it gives the proper asymp- 
totic scaling &t t ^ r^-] To avoid an explicit account 
of the detector quantum noise we will consider only pro- 
cesses at frequencies oj <C eV/fi; in particular, we assume 
max(nT-Me|,|iJ|)<eV^. 

Notice that due to electron charge discreteness and 
stochastic nature of tunneling, the total number n{t) 
of electrons passed through the detector is sometimes a 
more convenient magnitude to work with than the cur- 
rent I{t) = edn{t)/dt. In particular, we will use n{t) 
instead of I{t) for the Bayesian formalism derivation in 
the next section. 



III. DERIVATION OF THE BAYESIAN 
FORMALISM VIA "BLOCK" EQUATIONS 

The "conventional" ensemble-averaged equations for 
the qubit density matrix Pij{t), 

pn = -p22 = ~2 — Impi2, (7) 

£ H 

P12 = i ^/Oi2 + i-^ (pii - P22) - rdpi2, (8) 

do not take into account any information about the de- 
tector outcome and describe the effect of continuous mea- 
surement by the ensemble decoherence rate [13] 



2e 



(9) 



(Notice a relation Fdrm = 1/2; as will be seen later, this 
means that the detector is ideal). 



Equations (7)-(8) imply tracing over all detector de- 
grees of freedom, including the measurement outcome. 
An important step towards taking into account the mea- 
surement record was a derivation [13] of "Bloch" equa- 
tions for the density matrix plj{t) which is divided into 
components with different number n of electrons passed 
through the detector: 



Pn = -f + f Pn' - 2 I Imp?2 , 



P22 ~ P22 ~^ P22 
e e 

s H 
P12 = i ^ P12 + i ^ (Pii - 



2|^Imp5'2, 



P^22) 



h 



2e 



■P12 



+ ^P?2-^ 



(10) 
(11) 

(12) 



Eqs. (7)-(9) can be obtained from the Bloch equations 
using summation over n and relation pij = p^- . (Ab- 
sence of nondiagonal in n matrix elements is related 
to the assumption of large detector voltage [13].) 

Despite the Bloch equations carry the total number 
n of electrons passed through the detector, they can- 
not take into account the whole measurement record n(t) 
for a particular realization of measurement process. We 
should make a simple but important step for that: we 
should introduce a sufficiently frequent collapse of n(t) 
corresponding to a particular realization of the measure- 
ment record [4]. This idea is illustrated in Fig. 2(b). In- 
cluding "detector" into the quantum part of the setup, 
we anyway have to deal with a classical information, so 
we introduce a classical "pointer" which periodically (at 
times tk) forces the system "qubit+detector" to choose 
a definite value for n{tk). An obvious drawback of such 
construction is that it is absolutely not clear what should 
be a sequence of tk (in other words, how strongly the 
detector and pointer should be coupled). In general, 
the frequency of this collapse can depend on the phys- 
ical parameters of interaction between the measurement 
stage included in the "detector" Hamiltonian and the 
next stage. The only obvious fact is that if in an ex- 
periment we can record n{t) with some frequency, then 
the collapse should take place at least not less frequent. 
Still it is unclear how much more frequent. Fortunately, 
for a model described in the previous section the results 
do not depend on the choice of tk if Atk = tk — tk-i are 
sufficiently small, so the natural choice is Atk 0. 

Technically the procedure is the following. During the 
time between tk-i and tk the "qubit+detector" evolves 
according to the Bloch equations (10)-(12), while at time 
tk the number n is collapsed in the "orthodox" way [5]. 
This means that the probability P{n) of getting some 
n{tk) is equal to 



Pin) = p^,{tk) + p^2itk), 



(13) 



while after a particular Uk is picked, the density matrix 
pfj is immediately updated (collapsed): 
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PS(*fe + 0) = ^n,n. Pij{tk + 0) 
Pij{tk + 0) 



p'l^itk-o) + p^Utk-oy 



(14) 
(15) 



where (5„,„^ is the Kronecker symbol. After that the evo- 
lution is again described by Eqs. (10)-(12) with n shifted 
by rife until the next collapse occurs at i = i/t+i, and so 
on. This procedure is the m,ain step in the derivation of 
the Bayesian formalism. 

Let us discuss now the relation of this procedure to the 
classical Baycs theorem [10,11] which says that a posteri- 
ori probability P{Bi\A) of a hypothesis Bi after learning 
an information A {Bi form a complete set of mutually 
exclusive hypotheses) is equal to 



P{Bi\A) 



P{B,)P{A\Bi) 
Y.^P{Bk)P{A\Bk) 



(16) 



where P{Bi) is a priori probability of the hypothesis Bi 
(before learning information A) and P{A\Bi) is the con- 
ditional probability of the event A under hypothesis Bi. 

Assuming for a moment H = and £ = in the qubit 
Hamiltonian (so that the qubit evolution is due to mea- 
surement only), it is easy to find [4] that Eqs. (10)-(12) 
and our procedure (14)-(15) lead to the qubit evolution 



as 



Pii(tfe) = 



P22{tk) 



pii{tk-i)Pi{Ank) 



Pl2{tk) = Pl2{tk-l) 



pii(tfe_i)Pi(Anfe) + p22{tk-i)P2{Ank) 

P22{tk-l)P2{^nk) 

pii(tfe_i)Pi(Anfe) +p22(ife^i)P2(Anfc)' 

[pil{tk)p22{tk)V'^ 



(17) 



(18) 



(19) 



[pil{f^k^l)p22{tk-l)Y'^' 

Uk—rik-i is the number of electrons passed 



where An^ 
through the detector during time Atk and 

P,(n) = exp(-/,Aife/e) 



(20) 



is the classical Poisson distribution for this number as- 
suming either qubit state |1) or |2). One can see that the 
diagonal matrix elements pu exactly obey the classical 
Bayes formula (16), i.e. as if the qubit is really either in 
the state |1) or |2), but not in both simultaneously. Ac- 
tually, this fact is not much surprising because at least 
in some sense pu are the probabilities. Equation (19) is 
a little more surprising and says that the measurement 
preserves the degree of qubit purity ^12/(^11^22)^''^; for 
instance, a pure state remains pure during the whole mea- 
surement process. 

After introducing the main procedure (14)-(15), fur- 
ther derivation of the Bayesian formalism is pretty sim- 
ple and depends on whether we want to consider finite 
detector response, |A/| ~ /q = (/i -f /2)/2 or a weak 
response, |A/| <C lo- In the first case each event of tun- 
neling through the detector carries significant informa- 
tion and significantly affects the qubit state, so a rea- 
sonable "experimental" setup implies recording the time 



of each tunneling event. Then during the time peri- 
ods when no electrons are passing through the detector, 
the evolution is essentially described by the Bloch equa- 
tions (10)-(12) with n = 0, while the frequent collapses 
[Atk ^ inm{e/Ii,e/l2,h/H,h/e)] just restore the den- 
sity matrix normalization, leading to the continuous (but 
not unitary!) qubit evolution [4,9]: 



Pii = -P22 = -2 — lmpi2 P11P22, 

n e 



(21) 



ie iH AI 

P12 = ^ P12 + — [pn - P22) + —[Pn - P22)Pi2- (22) 



2e 



However, at moments when one electron passes through 
the detector, the qubit state changes abruptly (corre- 
sponding to Arife = 1 and Atk ^ in Eqs. (17)-(19)): 



pn{t + 0) 

p22{t + 0) 



hpuit-O) 



pl2{t + 0)=Pl2{t-0) 



hpll{t - 0) + l2P22it - 0) 

l-pii(i + 0), 

Pll{t + 0) P22{t 



0) 



Pii(t-O) P22(i-0) 



1/2 



(23) 
(24) 

(25) 



These abrupt changes arc usually called "quantum 
jumps" [9]. Notice that for Ji > I2 the jumps always 
shift the qubit state closer to |1) (because detector tun- 
neling is "more likely" for state |1)), while continuous 
nonunitary evolution shifts the state towards |2). On 
average the evolution is still given by conventional Eqs. 
(7)- (9). 

The case of a weak detector response, |A/| <C Iq, is 
more realistic from the experimental point of view. In 
this case the measurement time as well as the en- 
semble decoherence time are much longer than the 
average time e//o between electron passages in the de- 
tector. If the tunneling in the qubit is also sufficiently 
slow, h/H e/Io, we can completely disregard individ- 
ual events in the detector and consider the detector cur- 
rent I{t) as quasicontinuous. Then Eqs. (17) (18) for the 
evolution due to measurement only (neglecting H and e) 
transform into equations which again have clear Bayesian 
interpretation: 



Pll{t + T) 
P22{t + T) 



Pn{t)Pi{I) 



Pn{t)PiiI) + p22{t)P2{iy 

P22{t)P2(l) 
Pll{t)Pl{I) + P22{t)P2{iy 



where 



- 1 /■*+^ 

7 = - / I{t') dt' 



(26) 
(27) 

(28) 



is the detector current averaged over the time interval 
{t,t + T) and Pi{I) are its classical Gaussian probability 
distributions for two qubit states: 
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Piil) 



(27rD)i/2 



exp[ 



2D 



I D = So/2t, (29) 



(here 5*0 = 2eJo is the detector noise), while Eq. (19) 
essentially does not change. 

Differentiating Eqs. (26), (27), and (19) over time and 
including additional evolution due to H and e, we obtain 
the main equations of the Bayesian formalism: 



Pn 



-P22 = -2 ^ Im/>12 + Pll/922 [I{t) - lo], (30) 

ti bo 

.e .ff, 

Pl2 = 1 T Pl2 + 1 -r [Pll - P22 ) 



-{Pll - P22)^ [I{t) - lo] Pl2. 



(31) 



In each realization of measurement the noisy detector 
outcome I{t) is different; however, for each realization we 
can precisely monitor the evolution of the qubit density 
matrix, plugging experimental I{t) into Eqs. (30)-(31). 
Let us stress again that these equations show the absence 
of any qubit dccohcrcncc during the process of measure- 
ment. Pure initial state remains pure; moreover, initially 
mixed state gradually purifies during the measurement if 
H [4,8]. The gradual state purification is essentially 
due to acquiring more and more information about the 
qubit state from the measurement record. 

While the qubit state does not decohere in each in- 
dividual realization of the measurement, different mem- 
bers of the ensemble evolve differently because of random 
I{t). Averaging Eqs. (30)-(31) over random I{t) and us- 
ing the relation [which follows from Eqs. (13), (17)-(18), 
and (20)] 



m 



^0 = ^ (pii 



p22)+m, 



(32) 



where ^(t) is a zero-correlated ("white") random pro- 
cess with the same spectral density as the detector noise, 
~ So, wc obtain conventional Eqs. (7)-(8). There- 
fore, the ensemble- averaged decoherence in our model is 
just a consequence of averaging over the measurement 
outcome (similar conclusion is also valid in the finite re- 
sponse case). 

Notice that since I{t) contains the white noise con- 
tribution, Eqs. (30)-(31) arc nonlinear stochastic dif- 
ferential equations [23] and dealing with them requires 
a special care. The problem is that their analysis de- 
pends on the choice of the derivative definition. Two 
mainly used definitions are the symmetric derivative: 
p{t) = lim^^o[p(* + t/2) - p{t - t/2)]/t which leads to 
the so-called Stratonovich interpretation of the stochas- 
tic differential equations, and the forward derivative: 
p{t) = limr^o[p{t + t) — p{t)]/T {ltd interpretation). 
Usual calculus rules remain valid only in the Stratonovich 
form [23] , so the physical intuition works better when us- 
ing Stratonovich definition. However, Ito interpretation 



allows simple averaging over the noise and because of 
that is usually preferred by mathematicians. Since we 
derived Eqs. (30)-(31) by a simple first-order differentia- 
tion, we automatically obtained them in the Stratonovich 
form (keeping the second-order terms in the expansion, 
we can obtain different equations, depending on the def- 
inition of the derivative). Since sometimes Ito form is 
more preferable, let us translate Bayesian equations into 
Ito form using the following general rule [23]. For an 
arbitrary system of equations 



Xi{t)=Gi{^,t)+Fi{^,t)^{t) 



(33) 



in Stratonovich interpretation, the corresponding Ito 
equation which has the same solution is 



ii{t) = G.(x, t) + ^Yl + ^(*) ' 

k 

(34) 



where Xi{t) are the components of the vector x(t), Gt 
and Fi are arbitrary functions, and the constant is 
the spectral density of the white noise process ^{t). Ap- 
plying this transformation to Eqs. (30)-(31), we get the 
following equations in Ito interpretation: 



Pll 



-P22 

£ 



-2 — Im pi2 

n 

H 



2A/ , 
P11P22 i[t) , 



So 



P12 = 1 ^ P12 + i (pii 



P22) 



A/ (A/)2 

-(Pll - P22)-^P12 ^{t) Pl2 ■ 

bo ibo 



(35) 



(36) 



while the relation between pure noise ^ {t) and the current 
I{t) is still given by Eq. (32). Notice that the last term 
in Eq. (36) does not actually mean the single qubit deco- 
herence (pure state remains pure), but is just a feature 
of the Ito form [it directly corresponds to the ensemble 
decoherence after averaging over ^{t)]. 



IV. DERIVATION OF THE FORMALISM VIA 
CORRESPONDENCE PRINCIPLE 

Another derivation [8] of the Bayesian formalism for a 
single qubit can be based on the logical use of the cor- 
respondence principle [5], classical Bayes formula, and 
results of the conventional ensemble-averaged formalism. 
Even though this way lacks some advantages of the "mi- 
croscopic" derivation discussed in the previous section, it 
can be applied to a broader class of solid-state detectors, 
in particular, to the finite- transparency quantum point 
contact and (with some extension) to the single-electron 
transistor and SQUID. In this section wc will assume 
a double-dot qubit measured by a finite-transparency 
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QPC and treat the detector current I{t) as a quasicon- 
tinuous noisy signal that implies weak detector response, 
|A7| « Jo. 

Let lis start with a completely classical case when the 
electron is actually located in one of two dots and does 
not move, but we do not know in which one, so the mea- 
surement gradually reveals the actual electron location. 
This is a well studied problem of the probability theory. 
The measurement process can be described as an evolu- 
tion of probabilities (we still call them pn and P22) which 
reflect our knowledge about the system state. Then for 
arbitrary r (which can be comparable to Tm) the current 
/ averaged over time interval {t, t + t) [see Eq. (28)] has 
the probability distribution 

P(7) = Pll(i)Pl(7) + P22(t)P2(I), (37) 

where Pi arc given by Eq. (29) and depend on the detec- 
tor white noise spectral density which should not nec- 
essarily satisfy Schottky formula. After the measurement 
during time t the information about the system state has 
increased and the probabilities pn and P22 should be up- 
dated using the measurement result I and the Bayes for- 
mula (26) (27). This completely describes the classical 
measurement process. 

The next step in the derivation is an important as- 
sumption: in the quantum case with H — the evolu- 
tion of Pll and P22 is still given by Eqs. (26)-(27) because 
there is no possibility to distinguish between classical and 
quantum cases, performing only this kind of measure- 
ment. Even though this assumption is quite obvious, it 
is not derived formally but should rather be regarded as 
a consequence of the correspondence principle. 

The correspondence with classical measurement can- 
not describe the evolution of p\2', however, there is still 
an upper limit: \pi2\ < [piiP22]^^^ ■ Surprisingly, this in- 
equality is sufHcient for exact calculation of pi2{t) in the 
case of a QPC as a detector (we still assume H = 0). Av- 
eraging the inequality over all possible detector outputs 
/ using distribution (37) we get the inequality 

\{Pl2{t + T))\ < {\pi2{t + T)\) < {[p,,(^t + T)p22{t + r)]'/') 

= [pii{t)p22{t)]'^'exp[-{AlfT/4So] (38) 

[here the decaying exponent is a consequence of chang- 
ing Pll and P22 that reduces their average product]. On 
the other hand, from the conventional approach we know 
[16,18,24,25] that the ensemble-averaged qubit decoher- 
ence rate caused by a QPC is equal to = (AI)^ /ASq, 
where = 2e/o(l — T) and T is the QPC transparency. 
This means that inequality (38) actually reaches its up- 
per bound. This is possible only if in each realization of 
the measurement process an initially pure density matrix 
Pij{t) stays pure all the time, |pi2p = P11P22. Moreover, 
since the phase of pi2(i + t) should be the same for all 
realizations (to ensure |(pi2(i + T))| = (|pi2(t + t)|)), the 



only possibility in absence of a detector-induced shift of 

e is 

Pl2{t + T) ^ Pl2{t) i^^/R 

[pii(t + T)p22(i + r)]V2 [pn{t)p22{t)]y^ ^ ' 

(if the coupling with detector shifts e, we just have to use 
the shifted value). 

As the next step of the derivation, let us allow an ar- 
bitrary mixed initial state of the qubit (but still H = 0). 
It can always be represented as a mixture of two states 
with the same diagonal matrix elements, one of which 
is pure, while the other state does not have nondiago- 
nal matrix elements. Since nondiagonal matrix elements 
for the latter state cannot appear in the process of mea- 
surement and since the evolution of the diagonal matrix 
elements is equal for both states, one can conclude that 
Eq. (39) remains valid, i.e. for mixed states the degree of 
purity is conserved (gradual purification does not occur 
at = 0). The final step of the formalism derivation 
is differentiating Eqs. (26), (27), and (39) over time and 
adding (in a simple way) the evolution due to H. 

In this way we reproduce Eqs. (30)-(31). However, 
as seen from the derivation, now they are applicable to 
a broader class of detectors (which includes the finite- 
transparency QPC) for which = {AI)'^/4So. This 
relation can also be expressed as TdTm = 1/2 since Tm = 
25'o/(A/)^ for a weakly responding detector. As will be 
discussed in the next section, this is a condition of an 
ideal quantum detector. 

V. NONIDEAL DETECTORS 

The relation TjTm =1/2 which is valid for the mod- 
els of a tunnel junction and QPC as detectors, basically 
says that the larger output noise So of a detector leads 
to a smaller backaction characterized by ensemble deco- 
herence F^. This is quite expected from quantum me- 
chanical point of view (the faster we get information, the 
faster we should collapse the measured state). However, 
it is obviously not necessarily the case for an arbitrary 
solid-state detector; for example, the increase of output 
noise Sq can be due to later stages of signal amplifica- 
tion, which do not affect F^. In other words, it is easy 
to imagine a bad detector which produces a lot of both 
output and backaction noises. 

To take into account an extra detector noise, we can 
phenomenologically add a dephasing rate 7^ into the 
Bayesian equations: 

H 2 A/ 

Pll = -P22 = -2 -r IniPi2 + P11P22 [I{t) - h], (40) 
n bo 

.£ .H. 
P12 = 1 ^ P12 + 1 (Pii - P22) 

-(Pll - P22)^ [I{t) - lo] P12 - 7dPl2- (41) 
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This obviously increases the ensemble decoherence rate: 

(42) 



Td = —r^ h 7rf- 



4So 

A natural definition of a detector ideality (quantum effi- 
ciency) in this case is 



_ 1 7d 1 



(43) 



An upper limit for t] is 100% because of a fundamental 
limitation 



TdTm > 1/2, 



(44) 



which is a by-product of the Bayesian derivation for the 
case of quasicontinuous detector current and small differ- 
ence between noises Si and 5*2 (so that Si = S2 = So) 
- see inequality (38). [In the case of a detector with 
Si ^ 52 and possibility to observe each passing electron, 
Eq. (44) remains valid; however, a meaningful model 
would imply Poisson statistics (20) and therefore Schot- 
tky formula for the detector noise.] 

The extra dephasing 7^ in Eq. (41) can be interpreted 
[8] as an effect of extra environment or (mathemati- 
cally) as due to a second detector "in parallel" , the out- 
put of which is not read out (then we have to aver- 
age over possible outputs). It can be also interpreted 
as an effect of extra noise Sadd at the detector output, 
Sq = {AI)'^/4:Td + Sadd- In this latter case one can argue 
that the qubit is actually in a pure state and the evolu- 
tion of the diagonal matrix elements is actually different 
from what is given by Eq. (40), because the measured 
current I{t) is not the "actual" detector current. Yes, 
we would know the exact pure state if our amplifiers did 
not produce extra noise Sadd', however, since we do not 
have access to the "actual" detector current, we should 
perform averaging over the extra noise. It is easy to show 
that such averaged qubit density matrix (which is a den- 
sity matrix "for us") still satisfy Eqs. (40)-(41). 

Introduction of the detector ideality r] allows us to 
consider a continuous transition from the conventional 
ensemble-averaged result (7)-(8) to the "pure" Bayesian 
result (30)-(31). The effect of a pure environment can 
be considered as a measurement with an extremely bad 
detector, rj = 0. Technically it corresponds to AI = 
in Eqs. (40)-(41), transforming them into conventional 
equations. The case of a detector with very small effi- 
ciency, ry <C 1, can be treated in two steps: first, we an- 
alyze the effect of the decoherence term (7^ « Td), and 
then we use the classical (still Bayesian) analysis to relate 
the qubit density matrix and the measurement outcome. 
So, only for good detectors with the efficiency r] compa- 
rable to unity, the quantum Bayesian approach discussed 
in this paper is really necessary. Some people could argue 
that it is so difficult to create a solid-state detector with 



good quantum efficiency, that the Bayesian formalism is 
useless at the present-day level of technology. However, 
actually at present such detectors are becoming available. 
For example, the analysis of experimental data of the re- 
cent "which path" experiment [26] shows that their QPC 
had a quantum efficiency quite close to 100%. 

The account of the detector nonideality by introducing 
extra decoherence rate jd into Eq. (41) implicitly assumes 
the absence of a direct correlation between the output de- 
tector noise and the backaction noise affecting the qubit 
energy asymmetry e. However, such correlation is a typi- 
cal situation, for example, for a single-electron transistor 
as a detector [27]. In this case the knowledge of the noisy 
detector output I{t) gives some information about the 
probable backaction noise "trajectory" s{t) which can be 
used to improve our knowledge of the qubit state. Com- 
pensation for the most probable trajectory e{t) leads to 
the improved Bayesian equations for the SET in which 
Eq. (41) is replaced with 



n n 

AI 



P22) 



-{pii - 922)-^ [I{t) - h] pvi 
Jo 

-\-\K\I{t) - /o] P12 - 7dPi2, 



(45) 



where K = {de / dip)S / S^h, ip is the electric potential 
of the SET central electrode, and 5/^ is the mutual low- 
frequency spectral density between the current noise and 
ip noise. [Strictly speaking, Si^p in our notation is only 
the usual real part of the mutual spectral density, which 
reflects the detector "asymmetry", while the imaginary 
part can formally describe the detector response [16]. 
Also notice that since the small shift of the SET operat- 
ing point for two localized qubit states in general affects 
the energy £, it should be defined self-consistently in Eq. 
(45).] 

The dephasing rate 7^ should now satisfy equation 



7d = - 



{Alf K^Sq 
45o 4 



(46) 



to correspond to the the ensemble-averaged dynamics 
still described by Eqs. (7) -(8). 

The obvious inequality 7d > (in the opposite case 
the relation |y0i2p < P11P22 would be violated in a single 
realization of measurement) imposes a lower bound for 
the ensemble decoherence rate Fw: 



^ (Aiy K^So 
- 45o 4 



(47) 



which is stronger than the inequality TdTm > 1/2. 

Inequality (47) can be easily expressed in terms of the 
energy sensitivity of an SET. Let us define the output en- 
ergy sensitivity as e/ = {dl /dq)~^ Si/2C where C is the 
total SET island capacitance, dl/dq is the response to 
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the externally induced charge q, and we have changed the 
notation Sj = Sq for a more symmetric look of the for- 
mulas. Notice that e/ has the same dimension as U. Sim- 
ilarly, let us characterize the backaction noise intensity 
by = C Sip 1 2 and the correlation between two noises 
by the magnitude ei^ = {dI/dq)~^Si^/2. Since in ab- 
sence of other decoherence sources Yd = Sp{CAE/2eh)'^, 
where AE is the energy coupling between qubit and 
single-electron transistor [4,17], and using also the reci- 
procity property Ag' = CAE/e = ds/dip, we can rewrite 
Eq. (47) as 



e = (eje^ - ej^)'^^ > h/2. 



(48) 



This is a result known for 20 years [28] for SQUIDs (see 
also [16,29-33]). 

When the limit e = fi/2 is achieved, the decoherence 
rate 



7d = 



{All 
45/ 



- 1 



(49) 



in Eq. (45) vanishes, 7d = 0. In this sense the detector is 
ideal, fj = 1, where 



- _ 1 Id 



r{dl/dqf , [Sipf 



(50) 



even though it can be a nonideal detector (77 < 1) by the 
previous definition, = {dl /dqf / SjSp. 

Another possible definition of the detector efficiency in 
this case is 



V2 = 



h^idl/dqf 



c2 • 



Notice a simple relation. 



ri = r] = ri2 



2TdT„ 



(51) 



(52) 



in absence of correlation between noises of ^p{t) and I{t), 

Even though Eqs. (45)-(52) were derived for the SET 
as a detector, it is rather obvious that they are applica- 
ble to virtually any solid-state detector with continuous 
output (for a dc SQUID the current output should obvi- 
ously be replaced by the voltage output). In particular, 
the conclusion that reaching the quantum-limited total 
energy sensitivity e = h/2 is equivalent to the detector 
ideality, is quite general. 

As we already discussed, the tunnel junction and QPC 
at zero temperature (actually, for small temperatures 
-C eV) are theoretically ideal quantum detectors. 
The fact that a SQUID can reach the limit of an ideal de- 
tector follows from the results of Ref. [28] . A normal state 
SET is not a good quantum detector {rj <^ 1) at usual 
operating points above the Coulomb Blockade threshold 



[4,17]. However, its quantum efficiency improves when we 
go closer to the threshold [4,29] and becomes much bet- 
ter when the operating point is in the cotunneling range 
(below the threshold), in which case the limit of an ideal 
detector can be achieved [30,31]. Superconducting SET is 
generally better than normal SET as a quantum-limited 
detector and can approach 100% ideality in the super- 
current regime [32] as well as in the double Josephson- 
plus-quasiparticle regime [33] (possibly a threshold of a 
quasiparticle current is also a good operating point in 
this sense; however, this regime has been studied only 
for the current so far [34], but not for the noise). Finally, 
the resonant-tunneling SET [35] can reach ideality factor 
?72 = 3/4 at large bias and complete ideality, ?? = 772 = 1, 
in the small-bias limit. 



VI. SOME EXPERIMENTAL PREDICTIONS 
A. Direct experiments 

The Bayesian equations tell us that we can monitor the 
qubit evolution in a single realization of the measurement 
process using the record of the noisy detector output. In 
particular, for an ideal detector we can monitor the qubit 
wavefunction (except the overall phase) if the initial qubit 
state is pure or, for a mixed initial state, after monitoring 
for a sufficiently long time so that the gradual purification 
has enough time to produce a practically pure state. 

This prediction (and hence, the validity of the Bayesian 
equations) can in principle be tested experimentally. For 
example [8], let us first prepare the double-dot in the 
symmetric coherent state, pn = P22 = \pi2\ = 1/2, make 
H = (raise the barrier) , and begin measurement with a 
QPC [Fig. 1(b)]. According to our formalism, after some 
time T (the most interesting case is r r^) the qubit 
state remains pure but becomes asymmetric (pn ^ P22) 
and can be calculated with Eqs. (30) and (31). To prove 
this, an experimentalist can use the knowledge of the 
wavefunction to move the electron "coherently" into the 
first dot with 100% probabihty. (Notice that if the qubit 
is in a mixed state, no unitary transformation can end 
up in the state |1) with certainty.) For instance, experi- 
mentalist switches off the detector at t = t, reduces the 
barrier (to get finite H), and creates the energy difference 
e = [(1 -4|pi2(t)|2)1/2 _ l]i7Rcpi2(r)/|pi2(rj|2; then af- 
ter the time period At = [tt — arcsin(Impi2(T) hn/H)]/Q 
the "whole" electron will be moved into the first dot, that 
can be checked by the detector switched on again. [Here 
Q = (4iJ^ -|- e^y^^/h is the frequency of unperturbed 
coherent ("Rabi") oscillations of the qubit.] 

Another experimental idea [8] is to demonstrate the 
gradual purification of the double-dot density matrix. 
Let us start with a completely mixed (unknown) state 
(pii = P22 = 1/2, P12 = 0) of the double-dot qubit 
with finite H. Then using the detector output I{t) and 



8 



Eqs. (30)-(31) an experimentalist gradually gets more 
and more knowledge about the randomly evolving qubit 
state (gradual purification), eventually ending up with 
almost pure wavcfunction with precisely known phase of 
Rabi oscillations (we are not talking about the wavcfunc- 
tion phase, but about the phase of diagonal matrix ele- 
ments oscillations). The final check of the wavcfunction 
can be similar to that described in the previous para- 
graph. It can be even simpler, since with the knowledge 
of the phase of oscillations it is easy to stop the evolution 
by raising the barrier when the electron is in the first dot 
with certainty. Notice that if fast real-time calculations 
arc not available, the moment of raising the barrier can 
be random, while lucky cases can be selected later from 
the record of I{t). 

Direct experiments of this kind as well as experiments 
on quantum feedback control and on Bayesian measure- 
ment of entangled qubits (discussed in sections 7 and 8), 
are still too difficult for realization at the present-day 
level of technology. In the next two subsections we will 
discuss experiments which seem to be realizable (though 
very hard) at present. 

B. Spectral density of the detector current 

Naively thinking, a qubit with H ^ should perform 
coherent (Rabi) oscillations with frequency f2 and these 
oscillations should lead to an oscillating contribution of 
the detector current I{t). On the other hand, conven- 
tional Eqs. (7)-(8) seem naively to imply that the qubit 
eventually reaches a stationary state and no oscillations 
should be present in I{t) after a sufficiently long obser- 
vation. So, it is interesting to find what is the actual 
spectral density of the detector current Si{lo) [it is easier 
to measure this quantity experimentally, than to record 

m]. 

The Bayesian formalism predicts [21,19,22] the pres- 
ence of the spectral peak at the Rabi frequency fi, how- 
ever, the height of this peak cannot be larger than 4 times 
the noise pedestal. In particular for a symmetric qubit 

{8 = 0) 

So (w/fi)2 + (w2_f|2)2/f^2r2 • ^ ' 

Actually, an experimental confirmation of this formula 
would not be a direct verification of the Bayesian formal- 
ism, since Eq. (53) can be also obtained by other meth- 
ods, including the master equation method [25,21,16] and 
the method based on the Bloch equations [22]. 

C. Bell-type experiment 

Another experiment which also seems to be much eas- 
ier than the direct experiments but can unambiguously 



(a) 

on on 



del. A M qubit M det. B ^ 

Qa Qb 
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FIG. 3. (a) Schematic of the proposed Bell-type correlation 
experiment [36], in which a SCPB qubit is measured by two 
SETs during short time periods ta and tb shifted in time by 
r. The first measurement leads to an incomplete collapse of 
the qubit initial state (|1) -|- \2))/^/2 and affects the result of 
the second measurement, (b) The average result (Qb) of the 
second measurement for a selected result Qa of the first mea- 
surement. The sign and amplitude of Rabi oscillations depend 
on Qa, reflecting the change of the diagonal matrix elements 
of the qubit density matrix, (c) same as (b) if 7r/2 pulse is 
applied immediately after the first measurement. Now the 
phase of oscillations depends on Qa- The full-swing oscilla- 
tions (with amplitude of 0.5 in the ideal case) indicate a pure 
qubit state after the first measurement. 

test the; Bayesian formalism, is a Bell-type experiment 
in which one qubit is measured by two detectors [36]. 
An idea (Fig. 3) is to prepare the qubit in a coherent 
state (|1) + |2))/-\/2, then to switch on the first detector 
(A) for a relatively short time ta (so that the measure- 
ment is only partially completed), and to switch on the 
second detector (B) a little later. If the first measure- 
ment changes the qubit state according to the Bayesian 
formalism, then the second measurement can check this 
change. An output from a single run of the measure- 
ment are two charges Qa and Qb passed through two 
detectors. Performing the experiment many times and 
analyzing the correlation between Qa and Qb, one can 
recover the effect of the first measurement on the qubit 
state [36] (to check the change of the nondiagonal matrix 
element it is necessary to apply a 7r/2 pulse right after 
the first measurement) . The main advantage of this Bell- 
type experiment in comparison with the direct Bayesian 
experiments is that the wide bandwidth for the output 
signal is not necessary; instead, it is traded for the wide 
bandwidth of two input lines (switching detectors on and 
off), which is much easier to realize experimentally. 
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FIG. 4. (a) Schematic of the quantum feedback loop main- 
taining the Rabi oscillations of a qubit by synchronizing them 
with a classical harmonic signal, (b) Solid lines: the synchro- 
nization degree D as a function of the feedback factor F for 
several values of available bandwidth Tq"^. While synchroniza- 
tion can approach 100% for wide bandwidth, it worsens when 
Ta becomes comparable to the oscillation period T = 2n/Q. 
Dashed line (alsmost coinciding with the upper solid line): 
analytical result D — oxp(— C/32F). Dotted lino: synchro- 
nization degree for a direct feedback with Ta = T/10. (From 
Ref. [37].) 



VII. QUANTUM FEEDBACK CONTROL OF A 
QUBIT 

The Bayesian formalism can be used as a basis for the 
design and analysis of a quantum feedback control of a 
solid-state qubit. As an example, such feedback control 
can maintain for arbitrary long time the desired phase 
of a qubit Rabi oscillations, synchronizing them with a 
classical reference oscillator, even in presence of dephas- 
ing environment [4,37]. The overall idea is very close 
to a classical feedback loop [Fig. 4(a)]. The oscillating 
qubit evolution is monitored by a weakly coupled detec- 
tor (C = n{Alf/SoH < 1), the phase (j){t) of actual Rabi 
oscillations is compared with the desired phase (/'o(i)) and 
the difference signal Acf) is used to control the qubit bar- 
rier height. If qubit is slightly behind the desired phase, 
then H is decreased, so the oscillations run faster to catch 
up; if the qubit is ahead of proper phase, H is increased. 
It is natural to use a linear control: Hjb = iJ (1 — Fx Ai^), 
where F is a dimensionless feedback factor. 

The only difference of this loop from a classical feed- 
back is that even weakly coupled detector disturbs the 
qubit oscillations. However, the induced fluctuations of 
the oscillation phase are slow, and the information ob- 
tained from the detector happens to be enough to mon- 



itor the phase fluctuations and compensate them. The 
quantitative analysis [37] shows that in a limit of good 
synchronization and absence of extra environment the 
qubit correlation function Kz{t) = {z{t)z{t + r)) (here 
z = pu- P22) is given by 



COsflT 

Kz{t) = — - — exp 



16F 



-IFHtIU 



- 1 



(54) 



and does not decay to zero at t ^ 00. Correspondingly, 
the degree of the qubit synchronization, D = 2{ppd) — 1 
(here pd is the desired density matrix corresponding to 
ideal oscillations) is found iohe D = exp(— C/32F) and 
approaches 100% at F » C. 

The quality of the qubit oscillations synchronization 
decreases with the decrease of available feedback band- 
width [Fig. 4(b)]. It also decreases when the qubit 
is dephased by an extra environment. For a weak de- 
phasing rate 7e we found numerically [37] a dependence 
Dmax ^ 1 - 0.5de where 4 = 7e/[(A7)V45'o]. This 
means that the feedback loop can efficiently suppress 
the qubit dephasing due to coupling to the environment 
if this coupling is much weaker than the coupling to a 
nearly ideal detector. 

Besides the linear feedback Hfb = H{1 — F x At/)), we 
have also studied the "direct" feedback Hff,{t)/H — 1 = 
F{2[/(i) - Iq]/AI - cosflt} sinQt and found that it can 
also provide a good phase synchronization if F/C is close 
to 1/4 [Fig. 4(b)]. The direct feedback is much easier for 
an experimental realization because it does not require 
solving the Bayesian equations (40)-(41) in real time. 



VIII. BAYESIAN MEASUREMENT OF 
ENTANGLED QUBITS 

The Bayesian formalism has been generalized to a con- 
tinuous quantum measurement of entangled qubits in 
Ref. [38]. Suppose a detector is coupled to N entan- 
gled qubits. In the "measurement" basis there are 2^ 
states |i) of the qubits which correspond to up to 2^ dif- 
ferent dc current levels li of the detector (some of these 
currents can coincide, for example, if two or more qubits 
are coupled equally strong to the detector). It has been 
shown that the generalization of Eqs. (40)-(41) for this 
case is [38] 



Pij = y [Hqb, p]ij + Pij ^ ^ Pkk 



X{li-lk) + {l{t)- 



m 



lijpij, (55) 



where the first term describes the unitary evolution due 
to the Hamiltonian of qubits Hqb and 



(56) 
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while Eq. (32) is replaced by 



(a) 



(b) 



(57) 



Notice that there is no mutual decoherence (7^ = 0) 
between states \i) and |j) even for a nonideal detector 
if the corresponding classical currents coincide, 7, = Ij. 
This is because the detector noise cannot destroy the 
coherence between states which are equally coupled to 
the detector. 

Translating Eq. (55) from Stratonovich form into Ito 
form, we get 

Pij = y [Hgb, p]ij + pij (/(i) - ^ Pkklkj 

^ k 

X (/i + - 2 ^ pkkh) - (lij + ^^'Igf )pih (58) 



while in the ensemble-averaged equations the second 
term of Eq. (58) (which depends on I{t)) is averaged 
to zero. 

These Bayesian equations have been applied in Ref. 
[39] to the analysis of a simple setup (Fig. 5) in which 
a detector is equally coupled to two similar qubits (both 
qubits are symmetric, £a = £b = 0, and do not interact 
directly with each other). An interesting effect has been 
found in the case when the Rabi frequencies = 2Ha/h 
and fib = 2Hb/Ti exactly coincide. Then there are two 
possible scenarios of the two-qubit evolution due to mea- 
surement, starting from a general mixed state. Either 
qubits become fully entangled collapsing into the Bell 
state (Itoib)— |ioT6))/'\/2 (we call this process sponta- 
neous entanglement), or the state falls into the orthogo- 
nal subspace of the two-qubit Hilbert space. Experimen- 
tally these two scenarios can be distinguished by different 
spectral density Sii^jj) of the detector current [Fig. 5(c)]. 
In the case of Bell state, Si{ijj) is just the flat noise of 
the detector because the signals from two qubits compen- 
sate each other, while in the other scenario Si{u}) has a 
spectral peak at the Rabi frequency, which height is equal 
to 32?7/3 for a weakly coupled detector (C^ = C;, < 1) [39]. 

The probabilities of two scenarios are equal to the con- 
tributions of two subspaces in the initial state p(0); for 
the case of fully mixed initial state they are equal to 1/4 
and 3/4, respectively. The considered setup can obvi- 
ously be used for a preparation of the Bell state without 
direct interaction between two qubits. Notice that if the 
state collapsed into the orthogonal subspace, we can ap- 
ply some noise which affects and/or et and therefore 
mixes the two-qubit density matrix, and try the measure- 
ment again. In this way the probability 1 — (3/4)^ to 
obtain the Bell state can be made arbitrary close to 100% 
by allowing sufflciently large number M of attempts. 

In actual experiment the symmetry of the setup can- 
not be made exact. In this case the Bell state and the 
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FIG. 5. (a) Two qubits made of double quantum dots mea- 
sured by an equally coupled quantum point contact, (b) 
Similar setup made of single-Cooper-pair boxes measured by 

a single-electron transistor, (c) Two Monte Carlo realiza- 
tions of the two-qubit state evolution starting from the fully 
mixed state for a symmetric setup (pfi is the diagonal com- 
ponent of the two-qubit density matrix, corresponding to the 
Bell state {\]aib)- \ia]b))/-/2)- With probability 1/4 the 
qubits become fully entangled, pfi 1 ("spontaneous entan- 
glement"); then the detector output is a pure noise (upper 
inset). With probability 3/4 the state is gradually collapsed 
into the orthogonal subspace, pfi — » 0; then the detector sig- 
nal shows a spectral peak at the Rabi frequency Q, with the 
peak-to-pedestal ratio of 32/3. (Prom Ref. [39].) 



oscillating state are not infinitely stable and there will 
be switching between them. The calculations show [39] 
that the switching rate F b^o from the Bell state into the 
oscillating state is equal to Vb^o = iJ^^Y' /"^^d due to 
slightly different Rabi frequencies [F^ = 77^^ (A/)^/4S'o], 
Tb^O = {AC/C fTd/8 due to shghtly difl'erent coupling, 
and Tb^o = {"fa + "fb)/'2 due to an extra environment 
acting on two qubits separately. The rate of the return 
switching is 3 times smaller, Tq^b = Fb^o/3. Notice 
that in this case the averaged height of the Rabi spectral 
peak is equal to SrjSo, which is exactly twice as much as 
for a single qubit. 

Even though such experiment on spontaneous entan- 
glement is still extremely difficult for a realization, it 
should be noted that for the observation of the phe- 
nomenon the detector quantum efficiency r] should not 
necessarily be close to 100%; it should only be large 
enough to allow distinguishing the Rabi spectral peak 
with the peak-to-pedestal ratio of 3277/3. 
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IX. DISCUSSION 

In this paper we have reviewed the basic derivation and 
some apphcations of the Bayesian approach to continu- 
ous quantum measurement of sohd-state qubits. Even 
though this is a new subject for the soHd-state com- 
munity, many similar fornialisnis have been developed 
in other fields of quantum physics. Generally, this type 
of approach which takes into account the measurement 
outcome, is usually called selective or conditional quan- 
tum measurement. However, there is a rather broad va- 
riety of formalisms and their interpretations within the 
approach (for example, see reviews [40 43]). Some of 
key words related to this subject are: quantum trajec- 
tories, quantum state diffusion, quantum jumps, weak 
measurements, stochastic evolution of the wavefunction, 
stochastic Schrodinger equation, complex Hamiltonian, 
restricted path integral, quantum Bayes theorem, etc. 
The approach of conditional quantum measurements is 
relatively well developed in quantum optics. In particu- 
lar, the optical quantum feedback control has been well 
studied theoretically (see, e.g. [44 -48]) and was recently 
realized experimentally [49]. In relation to continuous 
quantum measurement of single systems, the quantum 
optics seems to be about 10 years ahead of the solid- 
state physics. However, the interest to this problem in 
the solid-state community has significantly increased af- 
ter the "which path" experiments [26,50]. Quite possibly 
it will be a rapidly growing field in the nearest future, 
especially because of its direct relation to the solid-state 
quantum computing. 

In this paper we have discussed two solid-state experi- 
ments which seem to be realizable (though very difficult) 
today. First, it would be interesting to measure the spec- 
tral density of the detector current when the measured 
qubit performs coherent oscillations, and compare exper- 
imental results with the theoretical prediction that the 
spectral peak in the best case is 4 times higher than the 
noise pedestal. Second, the Bell-type correlation exper- 
iment with one qubit measured by two detectors would 
be able to verify that the qubit state remains pure during 
the whole measurement process and show the possibility 
of monitoring the qubit evolution precisely. This would 
be the first step towards realization of the quantum feed- 
back control of solid-state qubits. A continuous monitor- 
ing of entangled qubits would be another very interesting 
experiment. Hopefully, the rapid progress of solid-state 
technology will make these experiments possible in a rea- 
sonably near future. 
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